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ABSTRACT 


Research  In  this  program  has  concentrated  on  the  generalized 
eigenvalue  problem  and  its  natural  extension  to  the  computation  of  the 
associated  canonical  form.  Furthermore,  there  has  been  an  extensive 
effort  to  study  the  matrix  equation  arising  in  control  engineering  such 
as  controllability  observability  decomposition  and  the  solution  of  the 
Riccati  equations.  In  particular,  error  bounds  for  the  computed  eigen¬ 
values  and  eigenvectors  of  the  generalized  eigenvalue  problem  have  been 
devised.  In  addition,  a  numerically  stable  algorithm  has  been  developed 
for  computing  the  orthonormal  bases  for  any  deflating  subspace  of  a 
regular  pencil.  A  method  has  been  developed  to  obtain  any  desired  order¬ 
ing  of  eigenvalues  in  the  quasi-triangular  forms. 


RESEARCH  OBJECTIVES  AND  STATUS  OF  THE  RESEARCH 


James  H.  Wilkinson 


I.  The  calculation  of  error  bounds  for  computed  eigenvalues  and 
eigenvectors 


In  practice  it  is  rare  for  the  complete  eigensystem  of  a  large  matrix 
to  be  required.  Commonly,  attention  is  focussed  on  a  few  key  eigenvalues 
and  eigenvectors.  An  algorithm  is  therefore  desirable  for  deriving  error 
bounds  for  a  single  eigenpair  (i.e.  value  and  vector)  without  requiring 
Information  on  the  remainder  of  the  system.  A  method  for  deriving 
rigorous  error  bounds  for  such  a  pair  X  and  x  was  developed  and  is 
about  to  be  published  in  Numexische  Mathematik  (a  preliminary  version  ap¬ 
peared  as  a  Stanford  report  [1]).  A  pleasing  feature  is  that  in  the  process 
of  deriving  the  bounds  an  Improved  eigenpair  was  derived  and  the  bound  was 
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for  this  eigenpair. 


A  less  pleasing  feature  was  that  the  bound  (and  indeed  the  true 
performance)  was  limited  by  the  condition  number  of  a  certain  matrix  C  . 
This  matrix  C  is  ill-conditioned  when  X  is  close  to  a  double  root. 
Multiple  roots  are  not  necessarily  ill-conditioned  and  hence  this  is  an 
inherent  weakness  of  the  method. 

This  weakness  has  been  overcome  by  determining  generators  of  the 
invariant  subspace  associated  with  clusters  of  eigenvalues.  This  material 
is  to  be  presented  as  in  invited  lecture  at  an  international  symposium 
at  ZUrich  in  honour  of  H.  Rutishauser  and  the  paper  [2]  will  be 
published  in  the  proceedings  of  that  symposium. 

Extensions  of  the  algorithm  which  make  it  computationally  more 
efficient  have  been  developed  by  J.  Dongarra  of  Argonne  National  Laboratory 
in  his  Ph.D.  thesis.  A  joint  paper  on  these  extensions  by  Dongarra,  Moler 
and  Wilkinson  is  in  preparation.  All  algorithms  apply  both  to  the  standard 
and  the  generalized  problem. 

II.  The  determination  of  the  distance  of  a  matrix  A  from  the  nearest 
defective  matrix  and  the  corresponding  problem  for  Ax  ■  Bx 

This  is  a  problem  of  considerable  interest  to  control  engineers. 

Early  work  in  connection  with  the  standard  problem  was  done  by  Kahan, 

Ruhe  and  Wilkinson.  The  techniques  developed  by  them  have  been  greatly 
improved  and  extended  to  the  generalized  problem  (the  more  important  case) . 
In  the  course  of  this  work,  rather  general  results  were  derived  which 
should  be  of  considerable  value  in  the  backward  error  analysis  of  eigen¬ 
value  algorithms.  This  work  is  to  be  presented  as  an  invited  lecture  at 
an  international  meeting  to  be  held  in  Winnipeg  in  October  and  the  paper 
[3]  will  be  published  in  the  proceedings. 
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III.  The  double  QR  algorithm  of  Frances 

This  algorithm  reduces  a  general  real  matrix  A  to  quasi- 
triangular  matrix  T  ,  which  has  two  by  two  matrices  on  the  diagonal 
corresponding  to  complex  conjugate  pairs  of  roots.  A  disadvantage  of 
the  algorithm  is  that  the  eigenvalues  are  not  ordered  in  any  systematic 
way  and  orderings  are  often  required  by  control  engineers. 

Algorithms  are  required  for  interchanging  the  order  of 

(a)  Adjacent  diagonal  elements 

(b)  A  two  by  two  block  and  an  adjacent  diagonal 

(c)  Two  adjacent  two  by  two  blocks. 

The  first  is  almost  trivial  (it  was  first  derived  by  Ruhe).  The  other 
problems  have  been  referred  to  briefly  by  Golub  and  Wilkinson  but  no 
solutions  were  given.  Methods  for  achieving  the  Interchanges  based  on 
orthogonal  similarities  have  been  developed  and  will  be  the  subject  of  a 
forthcoming  paper  when  they  have  been  fully  tested. 
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Paul  Van  Dooren 


I.  Computation  of  zeros  of  linear  multivariable  systems 

Several  algorithms  have  been  proposed  in  the  literature  for  the 
computation  of  the  zeros  of  a  linear  system  described  by  a  state-space 
model  {XI  -  A,B,C,D}.  Ue  have  devised  a  new  algorithm;  the  new  approach 
handles  both  nonsquare  and/or  degenerate  systems  without  difficulties 
whereas  earlier  methods  either  fail  or  require  special  tratment  for 
these  cases.  The  method  is  also  backward  stable  in  a  rigorous  sense. 

II.  A  generalized  eigenvalue  approach  for  solving  Riccatl  equations 

A  numerically  stable  algorithm  has  been  derived  to  compute  ortho- 
normal  bases  for  an  deflating  subspace  of  a  regular  pencil  XB-A  . 

The  method  is  based  on  an  update  of  the  QZ  algorithm,  in  order  to  obtain 
any  desired  ordering  of  eigenvalues  in  the  quasi-trlangular  forms 
constructed  by  this  algorithm. 

The  computation  of  deflating  subspaces  with  specified  spectrum  is  of 
crucial  importance  in  solving  Riccati  equations  arising  in  linear  system 
theory. 

Daniel  Boley 

I.  Computing  the  controllability-observability  decomposition  of  a  linear 
time-invariant  dynamic  system,  a  numerical  approach 

Various  numerical  properties  are  involved  in  computing  the  complete 
Controllability-Observability  (Kalman)  Decomposition  for  a  linear  time- 
invariant  dynamic  system,  of  the  form 

x  -  Ax  +  Bu 
y  ■  Cx 


-4- 


where  A,B,C  are  matrices,  and  u,  x,  y  are  vector  functions  of 

time.  In  particular,  we  are  investigating  the  numerical  stability,  the  cost 

and  the  particular  advantages  of  several  algorithms. 
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